High-temperature superfluidity of fermionic atoms in optical lattices 
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The experimental realizations of degenerate 
Bose |l| and Fermi atomic samples have 

stimulated a new wave of studies of quantum 
many-body systems in the dilute and weakly 
interacting regime. The intriguing prospective 
of extending these studies into the domain of 
strongly correlated phenomena is hindered by 
the apparent relative weakness of atomic inter- 
actions. The effects due to interactions can, 
however, be enhanced if the atoms are confined 
in optical potentials created by standing light 
waves The present Letter shows that these 

techniques, when applied to ensembles of cold 
fermionic atoms, can be used to dramatically in- 
crease the transition temperature to a superfluid 
state and thus make it readily observable un- 
der current experimental conditions. Depending 
upon carefully controlled parameters, a transition 
to a superfluid state of Cooper pairs antiferro- 



magnetic states [10] or more exotic d-wave pairing 
states [11] can be induced and probed. The re- 
sults of proposed experiments can provide a crit- 
ical insight into the origin of high-temperature 



superconductivity in cuprates [12|. 

An active search is now under way to implement a 
BCS transition of degenerate fermionic gases to a su- 
perfluid state analogous to superconductivity [ ji3fcf[4| . 
However, in free space or weakly confining atom traps 
the transition temperature to the superfluid state scales 
exponentially with interaction strength, ksT^ 8 « 
0.3££ oc exp[-7r/(2fc F |a s |)], with E f * cc the Fermi energy. 
For a dilute atomic gas the product of Fermi momen- 
tum and scattering length fcjr |a s | << 1, which makes the 
transition temperature exceedingly low. 

Atoms in potentials created by standing light waves 
(optical lattices) tend to localize near potential nodes 
thereby increasing the strength of effective interactions, 
whereas the kinetic energy provided by tunneling be- 
tween the different sites can be strongly suppressed [||. 
Very recently, fascinating experiments involving bosonic 
atoms in optical lattices [Q revealed a quantum phase 
transition from a superfluid to Mott insulating state 
Hp)] . Fermionic atoms confined in an optical lattice can 
undergo a phase transition to a superfluid state at a tem- 
perature that exceeds that of weakly confined atoms by 
several orders of magnitude. Attractive atomic interac- 
tions result in an s-wave pairing in which case fermionic 



atoms can undergo a BCS-type transition. In what per- 
haps is the even more intriguing prospective, fermionic 
atoms with repulsive interactions correspond to an ex- 
perimental realization of a Hubbard model that is widely 
discussed for strongly correlated electron systems such 
as high-T c cuprates [flGf . In particular, d-wave super- 
conducting states have been conjectured to exist in such 
systems, but so far this conjecture eluded rigorous con- 
firmation. We show that atomic systems with carefully 
controllable parameters and a variety of precise tools to 
detect the resulting phases can be used to provide a crit- 
ical insight into this outstanding problem. In essence, 
this approach can be viewed as an implementation of the 
pioneering ideas due to Feynman ]l7 [ for simulations of 
one quantum system by another. 

Consider an ensemble of fermionic atoms illuminated 
by several orthogonal, standing wave laser fields tuned 
far from atomic resonance. These fields produce a peri- 
odic potential for atomic motion in two (or three) dimen- 
sions of the form V(x) = Vq J2i=i cos 2 (kxi) with k the 
wave- vector of the light. The potential depth Vq is typ- 
ically expressed in the units of the atomic recoil energy 
Er = h 2 k 2 /2m. We will be interested in the situation 
in which there is roughly one atom per lattice site. Such 
atomic densities correspond to free-space Fermi energies 
on the order of Ep cc — (3/ir) 2 / 3 E R . At the same time the 
atoms can tunnel from one site to another. We assume 
that two kinds of atoms are present, differing by angular 
momentum or generalized spin (a = {T,4})- For suffi- 
ciently low temperatures the atoms will be confined to 
the lowest Bloch band, and the system can be described 
by a Hubbard Hamiltonian H|l| 



H 



C j,a C i,< 



where Ci ;CT are fermionic annihilation operators for lo- 
calized atom states of spin a on site i, m,„ — 
cta c i,&- We assume that the average occupation num- 
bers (ni,i;i) < 1. The parameter t corresponds to 
the tunneling matrix element between adjacent sites, 
t = Er(2JjJW)£? exp(— 2£ 2 ), and the parameter U = 
ERa s kyJ%/TT £ 3 characterizes the strength of the on-site 
interaction with £ = (Vo/Er) 1 / 4 . The sign of the scat- 
tering length a s determines the nature of atomic inter- 
actions: negative a s corresponds to attraction between 
atoms, whereas positive a s corresponds to repulsion. For 
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example, in the case of Li 6 ||,[|] both cases can be re- 
alized depending on the particular electronic states that 
are being trapped in a lattice. 

Consider first the situation corresponding to negative 
U. The effect of the lattice on the superfluid transition 
can be best understood starting from the limit of large 
tunneling i > \U\. Here, similar to the free-space tran- 
sition, the interaction strength is much weaker than the 
kinetic energy and the ground state of the system is then 
given by a "standard" BCS wave function, with an energy 
gap and transition temperature T c that depend upon t 
and u . For the lattice with filling fraction near unity 
((rif) + (ni) ~ 1) the Fermi energy is on the order of 
t, the density of states scales as N/t (with N being the 
total number of lattice sites) and the average two-atom 
interaction strength as U/N. Standard BCS theory can 
be applied to predict a critical temperature T c that for a 
3-D situation scales as fcsT c « 6texp(— 7t/|J7|) fL8| . An 
increase in the depth of the optical potential results in 
stronger atom localization and hence an increased inter- 
action strength U . At the same time, the tunneling t be- 
comes weaker. The combined effect of these two factors 
is a dramatic increase in T c due the exponential factor 
that competes with a very modest linear reduction due 
to a decreased Fermi energy. 

As the tunneling becomes comparable to the on-site 
interaction, the BCS picture is no longer valid. Due to 
strong attraction, atoms form pairs within single lattice 
sites. The entire system can then be considered as an 
ensemble of composite bosons. They can tunnel together 
at a rate ~ t /\U\, by virtual transitions via intermedi- 
ate singly occupied states. In this regime non-ordered 
pairs exist at high temperatures, whereas the superfluid 
state - a condensate of composite bosons - appears below 
k B T c ~ t 2 /\U\. Clearly, in this limit the increase in the 
potential depth will lead to a reduced mobility of pairs 
and hence a decrease in T c . The maximal critical tem- 
perature T c max is achieved at the crossover between the 
two regimes, when interaction and tunneling arc compa- 
rable (more precisely at U ~ 10i), which corresponds to 
a potential depth £ 2 s» l/21og[5v / 2/fc|a s |] and 

k B T™* » 0.3E$ cc k\a s \. (1) 

That is, the critical temperature for atomic fermions 
trapped in a lattice scales only linearly with the small 
parameter fc|a s |. This accurate result for the critical tem- 
perature is based on nonperturbative Monte-Carlo simu- 
lations of the fermionic Hubbard model . 

Several specific: approaches to achieve 1 the superfluid 
state can be considered. For example, the optical po- 
tential can be adiabatically turned on, starting from a 
weakly confined Fermi-degenerate mixture of the two 
atomic states of appropriately chosen density. In this 
procedure the atomic quasi-momentum is approximately 
conserved but the band-structure associated with the pe- 
riodic potential changes, resulting in a non-equilibrium 



distribution, with an effective temperature Tf different 
from the initial Tj n . The final temperature Tt can eas- 
ily be estimated from the relation X)k e k/( e k/^i«) = 
Ek £ k/(ek/T/), where f(x) = l/(e x + 1) is the Fermi- 
Dirac distribution function, e£ = fc 2 /2m — /i° is the 
original dispersion of atoms in free space, and et = 
— 2t(cos(A; x a) + cos(k y a) + cos(fc z a)) — \i is the dispersion 
in a tight-binding model with a = n/X the lattice period. 
Two important processes determine Tf. First of all, the 
presence of the lattice makes the system anisotropic and 
hence changes the shape of the Fermi surface. This re- 
sults in an effective heating of the system as Vo increases 
from zero to about E^. As Vq is increased even further 
the shape of the Fermi surface remains approximately 
the same and only the Fermi velocity (or effective atomic 
mass) changes leading to a reduction of the effective tem- 
perature within the lowest Bloch band (see Fig. 1). This 
suggests that it is optimal to turn on a weak lattice poten- 
tial while the fermionic sample is in contact with a cool- 
ing reservoir (e.g. atomic BEC). Although the resulting 
phase space density will not be significantly altered, this 
loading procedure will allow to avoid the heating which 
is present at the initial stages of creating an optical lat- 
tice. When Vo ~ En (point C in Fig. la) the system is 
decoupled from the reservoir and the lattice potential is 
increased until the transition to the superfluid phase is 
reached. 

To control precisely the resulting quantum phase (es- 
pecially in the situations described below) an accurate 
manipulation of the filling fraction may be important. 
This can be achieved, for example, if atoms with three 
internal states are used. If a dense, degenerate ensemble 
is prepared in a state that is not affected by the opti- 
cal lattice, a laser driven Raman transition into a pair of 
trapped spin states can be used to produce exactly one 
atom (or its fraction) per each lattice site. The essential 
idea of this approach is that the energy shifts associ- 
ated with atom interactions and the Pauli principle can 
be used to block the transitions into states with more 
than one atom per lattice site. In this case an effective 
^7r-pulse will result in a filling fraction of v with an un- 
certainty that scales at most with the inverse size of the 
lattice. 

Let us now consider the implications of these results in 
the light of present experimental possibilities. The rele- 
vant temperature calculated numerically for Li 6 and 
K 40 M atoms is shown in Fig. 1. In this figure we con- 
sider a Li 6 atomic sample of a very modest density corre- 
sponding to a unity filling in an optical lattice produced 
by CO2 laser (A — 10/im). For such densities, a dramatic 
increase in the critical temperature is possible. Note in 
particular that a phase transition can be achieved start- 
ing from an initial temperature of about 0.1E^- rcc . The 
CO2 lattice has the additional advantage of exceptionally 
long lifetimes, which should be sufficient to achieve the 
transition even for relatively low energy scales involved. 
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Another scenario is to trap Li atoms in an optical lattice 
created by a Nd:YAG laser A ~ 1.06/im. Although in this 
case the densities of 10 12 — 10 13 cm -3 will correspond to a 
filling fraction slightly less than unity, the resulting criti- 
cal temperature can still be in the range of 0.1E F CC . The 
inset shows a diagram for K 40 atoms trapped in a similar 
lattice. As indicated by the two cases presented in Fig. 1 
the transition to a superffuid state is expected to occur 
for the same initial temperature if cooling due to adia- 
batic switching on of the lattice is taken into account. 
Therefore, the value of the maximal initial temperature 
at which a phase transition can occur is almost indepen- 
dent of the scattering length and corresponds to about 
one tenth of the free space Fermi energy. 

Before proceeding we point out that in contrast to the 
approaches that are based on increasing the scattering 
length or atom densities, which result in a very inter- 
esting regime of BCS-BEC crossover [^9|-^T|, an optical 
lattice does not lead to an enhancement of inelastic loss 
processes but rather their suppression, as no more than 
two atoms can ever occupy the same lattice site. We also 
note that this implies an extremely low spin flip rate be- 
tween different atomic states, thereby allowing compara- 
tively large magnetic fields to tune the scattering lengths 
without significantly perturbing the relative fractions of 
n± and n-\. It is also important to point out that al- 
though the free-space ensembles at densities n|a s | 3 > 1 
can no longer be considered as weakly interacting gas 
p9| , p0[ , the Hubbard model nevertheless remains a valid 
description for the lattice case even in the regime of very 
strong confinment. Here corrections to the on-site inter- 
action can be accurately derived from known molecular 
potentials. Finally we emphasize again, that in contrast 
to a weakly confined Fermi-gas, the critical temperature 
for the optical lattice filled with attractive atoms (Fig. 1) 
can be predicted very accurately even in the most inter- 
esting, intermediate, regime t ~ \U\, since the behaviour 
of the Hubbard model for this case is by now very well 
understood. 

It is intriguing to consider possible extensions of the 
above ideas to a situation in which different atoms repel 
each other (o s > 0). This is realized for | J,) = \F — 
9/2, m F = 9/2) and | f) = \F = 9/2, m F = 7/2) states 
of atomic K at zero magnetic field U. In this case it is 
energetically unfavorable for two atoms to be on the same 
lattice site. However, adjacent atoms can virtually tun- 
nel to the same site. This process lowers the total energy 
of two atoms in adjacent wells, thereby creating effective 
hard-core attractive interactions of different spins. When 
the filling fraction of the lattice is close to one, this leads 
to a ground state in which adjacent sites are always occu- 
pied by atoms with alternating spins (see Fig. 2), i.e. an 
antiferromagnetic phase. 

What happens when the filling fraction is smaller than 
one? This question lies at the heart of the present de- 
bates on the nature of high-temperature superconducting 



cuprates. Clearly, a strong on-site repulsion makes it un- 
favorable for atoms to bind into the usual s-wave pairs, in 
which the probability to be at the center of orbit should 
be maximal. However, it has been conjectured [^l) that 
anisotropic d-wave pairs can be formed, which can result 
in a d-wave supcrfluid phase capable of explaining many 
of the observed properties in cuprates. 

d-wave Cooper pairs can be thought of as spin singlet 
angular momentum I = 2 Cooper pairs. In the absence 
of the lattice potential all five l z — ±2, ±1, components 
are degenerate. A cubic potential splits these states ac- 
cording to the representations of the cubic point group, 
and we find triply and doubly degenerate Cooper pairs 
. When the symmetry of the crystal is lowered even 
more, as e.g. in a two dimensional lattice, the degenera- 
cies between various Copper pairs are removed further, 
and in the x,y directions the lowest energy Cooper pair 
is described by the momentum space wave function |24[] 
A(p) = cos(p x a) — cos(p y a). Although so far the exis- 
tence of d-wave superconductivity in the repulsive Hub- 
bard model eluded rigorous confirmation, we believe that 
these ideas can now be tested experimentally in ensem- 
bles of fermionic atoms. For example, Fig. 2 shows a 
phase diagram for the system of repulsive atoms in two 
dimensions calculated using the FLEX approxima- 
tion [p5[ . Although the resulting T c is believed to be 
somewhat lower than in the s-wave case, this calculation 
suggests the existence of the d-wave phase for feasible 
atomic temperatures and densities. 

We next consider several approaches that can be used 
to detect and accurately probe the resulting quantum 
phases. For example, interference of the atoms released 
from the lattice has been used to probe the superfluidity 
of bosons [Q. Due to the exclusion principle, identical 
fermions can never be in the same quantum state. This 
implies that in the degenerate regime atomic interference 
patterns due to each momentum state will be superim- 
posed, resulting in real space interference peaks which 
reflect the shape of the Fermi surface. With the appear- 
ance of pairing, the atomic momentum in the pairs be- 
comes on the order of Planck's constant divided by the 
size of the pair. As a result, the momentum distribution 
is no longer a sharp step function, which will be directly 
reflected in the interference pattern as shown in Fig. 3. 

In order to directly detect superfluidity of the pairs, 
photoassociation spectroscopy can be used. Weakly 
bound Cooper pairs can be converted into molecules by 
using a laser-induced transition into a bound molecular 
state. The interference pattern of the released bosonic 
molecules will then provide extremely sharp peaks due 
to the presence of a supcrfluid fraction, in direct analogy 
to ref. §. 

The spectrum of elementary excitations also provides 
an accurate probe for the nature of the quantum phase. 
It can be measured in a system of cold atoms by excit- 
ing the motional states of atoms using laser pulses. For 
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example, atoms can experience Bragg scattering off two 
non-collinear laser beams, provided that the frequency 
difference 8ui of the lasers matches the resonance fre- 
quency of elementary excitation with momentum q deter- 
mined by the angle between two lasers (Fig. 4) p(| . This 
technique provides a direct measurement of the density- 
density correlation function. By monitoring the num- 
ber of Bragg-scattered atoms as a function of Sui, the 
presence of a superfluid phase can be detected unam- 
biguously: One observes a sharp peak due to a collective 
Bogoliubov mode, which is separated from a broad fea- 
ture (corresponding to quasi-particle excitations) by an 
energy gap (see Fig. 4). In the normal state, on the other 
hand, the continuum excitations are gapless and the col- 
lective mode is not visible due to strong damping. 

With this technique, it is also possible to detect d- 
wave superfluidity. For example, Fig. 5 shows the onset 
frequency of the quasiparticle continuum corresponding 
to the d-wave superfluid phase. A strong anisotropy, 
together with a vanishing gap for certain momenta can 
provide unambiguous evidence for the presence of such a 
phase. 

It should be noted that in addition to the effects de- 
scribed here, a number of other intriguing avenues can 
be considered. Atomic systems make it possible to ac- 
curately probe the effects of dimensionality. Varying 
for example the depth of the trapping potential in one 
direction while keeping the other parameters fixed, it 
should be possible to observe a transition between a two- 
dimensional and a three-dimensional situation. One can 
also envision a situation in which more than two internal 
atomic states are trapped in the optical lattice. Such a 
system is expected to result in a variety of new quantum 
phases such as superconductivity based on quadruplets 
rather than pairs. Likewise, these systems might allow 
for accurate studies of the effect of dissipation and de- 
coherence on macroscopic quantum phases. Another in- 
teresting phenomenon that may be studied by creating 
a system with unequal densities for the two spin states 
is the possibility of the Fulde-Ferrel-Larkin-Ovchinnkov 
(FFLO) state [|27|]28[ with a modulated superfluid order 
parameter. 
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FIG. 1. Left: Attractive fermionic atoms in optical lat- 
tices can undergo pairing resulting in a BCS transition into 
a superfluid state. Right: Critical temperature for a tran- 
sition of Li 6 atoms (circles) into the superfluid state as a 
function of the optical lattice depth in a three-dimensional 
CO2 lattice. Li atoms in the states with different orienta- 
tion of nuclear spin | 1) = \F = 1/2, = —1/2) and 

f) = \F = 1/2, mF = 1/2) are considered at a magnetic 
field of ~ 0.1 T which corresponds to a scattering length of 
a s ~ —2.5 x 10 3 ao in atomic units. The absolute energy 
scale is given by t/% ~ 0.5kHz at the phase transition. For 
the same values of the scattering lengths and densities the 
free-space BCS formula gives for the temperature of the su- 
perfluid transition T c ° = 1.6 x lQ- 12 E% cc /k B for Li 6 . The 
inset shows the analogous plot for K 40 atoms in a Nd:YAG 
lattice at half filling. For K 40 we took two states characterized 
by total angular momentum \ 1} = \F = 9/2,mF = —9/2) 
and I f) = \F — 9/2, tuf ~ —7/2) at a magnetic field above 
a Feshbach resonance and choose a s ~ —2. x 10 2 ao. As dis- 
cussed in the text, turning on an optical lattice also changes 
the temperature of the system. It is most advantageous to 
turn on a weak lattice potential Vq/Er ~ 1 while the atoms 
are cooled. The cooling is then switched off (point C in the 
figure) and the lattice depth is increased adiabatically, reduc- 
ing the effective temperature (shown by the dashed line for 
an initial temperature Tin = Q.Q7E l p ea /ks)- 




FIG. 2. In the case of repulsion it is energetically unfavor- 
able for two atoms to be the on the same lattice site. However, 
virtual tunneling creates effective hard-core attractive inter- 
actions of different spins on nearby sites. This may result 
either in antiferromagnetic (upper left) or d-wave superfluid 
phases, depending upon the filling fraction n. In three di- 
mensions and at filling fraction n = 1 the maximum critical 
temperature for antiferromagnetic ordering is identical to the 
superfluid case, and of the order 0.1E t J :cc . The d-wave sym- 
metry of the superfluid is related to the "parent" insulating 
phase (see lower left). The phase diagram for repelling Li 
atoms in a 2D lattice, calculated using the FLEX approxima- 
tion, is shown on the right hand side. Adiabatic cooling due 
to switching on of the lattice has been taken into account. 
The situation we have in mind corresponds to a very large 
potential depth along one direction and an equal, finite depth 
along the other two directions. This results in a set of weakly 
coupled 2D lattices. Within each lattice tunneling is given by 
t. We note that the repulsive Hubbard model may also have 
phase separation for some filling factors, which corresponds 
to immiscibility of the two spin species of the atoms. 




FIG. 3. Central part of the atomic interference pattern af- 
ter a free expansion for t = 500 x (Tl/Er). We are considering 
a 10x10 optical lattice (a denotes the lattice constant). Left: 
normal state at half filling. The sharp edges of the interfer- 
ence peaks reflect the atomic momentum distribution: in the 
normal state the peaks will become sharper with decreasing 
temperature - a direct consequence of Fermi-statistics. Right: 
BCS state at half filling, with a gap A/t — 0.6 corresponding 
to an interaction strength U jt w —2.5. 
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FIG. 4. Left: Bragg scattering of atoms off laser beams 
with frequency difference 8lj and wave-vector difference 
<lx — q y — O.Itt for attractive fermions at filling n = 0.6 
and U/t — —2.5, corresponding to an s-wave gap A/t m 0.45. 
The spectrum consists of a sharp collective mode and a contin- 
uum at higher frequencies. The sharp resonance corresponds 
to long-wave coherent excitations of the condensate. The 
continuum arises due to the breaking of Cooper pairs, which 
results in quasiparticle excitations. In the inset we show the 
dispersion of the collective mode (solid line) and the onset 
frequency of the continuum (dashed line). Right: schematic 
picture of the two-photon process involved. Qi and SI2 are 
the Rabi rates of the optical fields. 
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FIG. 5. Probing d-wave pairing at filling n < 1 via Bragg 
scattering. Left: schematic diagram of the Fermi surface 
(solid line) and the g-dependence of the Cooper pair wave 
function A(q x ,q y ) (dashed line). At the four nodal points 
shown by black dots, the wave function and the quasiparti- 
cle excitation energy vanishes. For the special wave vectors 
connecting these points, the density response is gapless (black 
spots in the right figure). Right: onset frequency LOmi n (q x , q y ) 
of the quasiparticle continuum, dark regions corresponding to 
low frequencies (vanishing gap). 
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